PROCEEDINGS 

OF SCIENCE 



AGN emission processes of NGC 4945 in the X-rays 
and y-rays 



Marie-Luise Menzel j Volker Beckmann, and Fabio Mattana 

Frangois Arago Centre, APC, Universite Paris Diderot, CNRS/IN2P3, CEA/Irfu, Observatoire 
de Paris, Sorbonne Paris Cite, 10 rue A. Domon et L. Duquet, 75205 Paris Cedex 13, France 
E-mail: }nl . men z el @t -online . de 



NGC 4945 has an outstanding role among the Seyfert 2 active galatic nuclei (AGN) because it 
is one of the few non-blazars which have been detected in the y-rays. Here, we analyse the high 
energy spectrum using Suzaku, INTEGRAL and Fermi data. We reconstruct the spectral energy 
distribution in the soft X-ray to y-ray domain in order to provide a better understanding of the 
processes in the AGN. We present two models to fit the high-energy data. 

The first model assumes that the y-ray emission originates from one single non-thermal com- 
ponent, e.g. a shock-induced pion decay caused by the starburst processes in the host galaxy, 
or by interaction with cosmic rays. The second model describes the high-energy spectrum 
by two independent components: a thermal inverse Compton process of photons in the non- 
beamed AGN and a non-thermal emission of the y-rays. These components are represented 
by an absorbed cut-off power law for the thermal component in the X-ray energy range and 
a simple power law for the non-thermal component in the y-rays. For the thermal process, 
we obtain a photon index of V = 1.6, a cut-off energy of £ C ut ~ 150keV and a hydrogen 
column density of N^ d ~ 6 x 10 24 cm~ 2 . The non-thermal process has a photon index of 
r = 2.0 and a flux of fb.i-iooGeV = 1-4 x 10~ n ergcm~ 2 s~'. The spectral energy distribu- 
tion gives a total unabsorbed flux of f2keV-i00GeV ~ 5 x 10 _10 ergcm _2 s _1 and a luminosity 
of £2keV-i00GeV ~ 9 x 10 41 ergs -1 at a distance of 3.7 Mpc. It appears more reasonable that the 
y-ray emission is independent from the AGN and could be caused e.g. by shock processes in the 
starburst regions of the host galaxy. 
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1. Introduction 

NGC 4945 is a nearby galaxy (z = 0.0019 or 3.7 Mpc) in the southern hemisphere at RA = 
13 h 05 m 27 s and DEC = -49°28 m 06 s (J2000.0) which is expected to belong to the Centaurus 
galaxy group [[Tj]. It is an almost edge-on spiral galaxy (inclination angle i = 78° [§]) of type 
SBcd or SABcd. 

NGC 4945 exhibits a prominent dust lane crossing its plane and hosts a bright and compact 
core. It is an example of a galaxy with a composite nature of an Seyfert 2 active galactic nuclei 
(AGN) and a star-forming region [^]. It shows AGN X-ray emission as well as a strong infrared 
component originating from starburst activity. NGC 4945 is the brightest AGN above 100 keV 
and one of the best studied Compton-thick sources with an intrinsic column density of the order of 
iV H = 5 x 10 24 cm~ 2 [Q]. From the definition of radio loudness as R* = f 5G uz //optical NGC 4945 
classifies as a radio quiet object, with an optical core brightness of ~ 9.5 mag and a radio flux of 
~ 3 Jy, resulting in R* ~ 4.5. Recently, Teng et al. ^ applied a different measure for the radio 
loudness, based on the radio to X-ray ratio. Here, NGC 4945 would qualify as slightly radio loud 
when compared to other Seyfert galaxies (by using the X-ray or hard X-ray flux). 

At high energies NGC 4945 is one of the few non-blazars which have been detected in the y- 
rays. This detection came as a surprise, because non-beamed AGN and starburst galaxies were not 
expected to produce significant emission at these high energies. This also motivated our study to 
investigate the spectral energy distribution of NGC 4945 in greater detail. Here we use unpublished 
hard X-ray data of NGC 4945 from the INTEGRAL and Suzaku X-ray satellites, and published 
results from Fermi/LAI. We analyze the combined X-ray and gamma-ray spectra in order to test 
different models to explain the processes driving the overall emission. In considering the star 
formation in the host galaxy and the processes in the AGN core, the physical connection between 
the emission processes of the X-ray to the y-ray spectrum is investigated. On one hand, thermal 
inverse Compton is assumed to dominate the spectrum of Seyfert galaxies in the X-ray domain, 
where the contribution of the starburst is expected to be faint. On the other hand, both the starbust 
activity (e.g., via collisionless shocks) and the AGN can contribute to the y-ray yield. It is important 
to disentangle these components in order to clarify the physical mechanisms dominating in AGN. 

2. Data analysis 

The spectrum of NGC 4945 which was analyzed in this work has been observed by instruments 
on-board Suzaku, INTEGRAL and Fermi. In Table [T], we give the characteristics of the observations. 

From Suzaku, we used data of XIS (X-ray Imaging Spectrometer) and HXD-PIN (Hard X- 
ray Detector) based on the observation of 15 January 2006. The XIS total spectrum has been 
derived by summing up the XIS-0, XIS-2 and XIS-3 spectra. The background-corrected spectra 
have been extracted from a 90" radius region centered on the source and data in the 4.5 - 8.5 keV 
band have been used. The spectrum from the non-imaging HXD detector have been background 
and dead time corrected following the standard analysis. The cosmic X-ray background (CXB) 
is approximately 5% of the NXB (non X-ray background) component and was estimated with a 
high-energy cut-off power law by using the HXD-PIN response for diffuse emission in XSPEC. 
HXD-PIN data have been extracted in the 15-40 keV band and the spectrum has been rebinned 
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mission 


instrument 


net exposure time 


energy range used 


Suzaku 


XIS (XO, X2, X3) 


285.3 ks 


4.5 - 8.5 keV 


Suzaku 


HXD-PIN 


84.9 ks 


15 -40keV 


INTEGRAL 


JEM-X 


152.5 ks 


10 - 35 keV 


INTEGRAL 


ISGRI 


67.2 ks 


40 - 250 keV 


INTEGRAL 


SPI 


530 ks 


25 - 70 keV 


Fermi 


LAT 


36.8 Ms 


150 MeV - 2.5 GeV 



Table 1: Observation log of the data used in this study. 



to have at least 300 counts per channel. In addition, we used data from the instruments SPI, ISGRI 
and JEM-X onboard the INTEGRAL satellite. Here we summed up the available data since the start 
of the mission in late 2002. A significant spectrum was extracted from SPI data in the 25 - 70 keV 
band. IBIS/ISGRI data showed some calibration problems at low energies and thus we used data 
only within 40 - 350 keV, and JEM-X spectra have been extracted from 10 - 35 keV. 

We used Fermi Large Area Telescope (LAT) observations covering August 2008 until end of 
March 201 1. The data have been extracted from a 3.3° location around the AGN in order to avoid 
contamination by other sources in the region. The dead time corrected data have been rebinned 
to 5 energy bins in the range 150 MeV - 2.5 GeV. Because of the low-statistics in source counts, 
we applied a simple power law in the gtlike analysis of the Fermi/LAI data. The background 
has been modeled by an extragalactic component with a flux of F extraga i = 3.1 x 10 _4 phcm _2 s _1 
and a Galactic background of F Gal = 4.8 x 10" 4 phones" 1 . As result, we obtain a photon in- 
dex of T = 2.3 ±0.1, a revalue of 92.9 (about 9.6a) and a flux of F (>100 MeV) = (2.2 ±0.5) x 
10~ 8 phcm~ 2 s _1 . Our fit parameters are consistent with the results derived by Lenain et al. [0]. 

3. Spectral energy distribution 

In order to interpret the multiwavelength spectrum in the 10 keV - 2.5 GeV range we focus 
on two models to explain the X-ray to y-ray spectrum: a) a non-thermal model or b) a combined 
thermal and non-thermal model. The analysis was performed using the XSPEC spectral fitting 
package version 12.5.1. 

First, we will assume that the X-ray to y-ray emission originates from one single non-thermal 
component. This component could be explained e.g. by mutiple shocks caused by the starburst 
processes in the host galaxy, or by interaction with cosmic rays. A simple absorbed power law 
model (wabs*powerlaw) gives a photon index of T = 2.1 ±0.1, a hydrogen column density 
of 7V H = (6.8 ±0.6) x 10 24 cm~ 2 with a goodness of fit of Xv = 1-06 for 66 degrees of freedom 
(d.o.f.). Adding an exponential cut-off to the data (wabs*cutof fpl) gives a photon index of 
r = 1.5 ±0.4, with an exponential cut-off of E cut < 478 keV. This model is a poorer representation 
of the data (with # 2 = 1.2 for 65 d.o.f.), and in addition it is not able to reproduce the LAT data 
points. 

The high-energy spectrum can also be described by two independent components. The first 
component is assumed to result from the photons from the thermal inverse Compton processes in 
the non-beamed AGN. The X-ray component is fitted by a cut-off power law (wabs*cutof fpl) 
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parameter 


value 




1.6 ±0.4 


(10 cm z ) 


7 (fixed) 


kT (keV) 


3+ 3 


E cut (keV) 


150 (fixed) 


N^ d (10 23 cm~ 2 ) 




r 2 


2.0+- 


%l (d.o.f.) 


1.16(98) 



line 


energy [keV] 


EW [eV] 


Fe I Kai 
Fe I Kg!2 
Fe xxv Ka 
FeiK/3 
NiiKa 


6.40 ±0.01 

6 - 25 -0.06 

6.68 ±0.02 

7 04 + 02 
'• u -0.03 

7.47 (fix) 


400+- 
80^ 
1403 
130 ±30 
70 ±50 



Table 2: Parameters of the fit. Errors have been determined on the X-ray spectrum without the Fermi/LKT 
data. For the iron line measurements, only the Suzaku/XIS data have been used. 



described by the photon index Ti, the cut-off energy E cut and the hydrogen column density Nu- 
The second component, responsible for the emission at y-rays is non-thermal and thus fit by an 
additional simple power law described with photon index F2- The fit results show that wabs 
*cutof f pl+powerlaw is the best model for the X-ray spectrum {% 2 = 1.01, 63 d.o.f.). As- 
suming an absorbed cut-off power law leads to a photon index of Y\ = 1 .6 ± 0.4, an upper limit for 
the cut-off energy of E cut < 230 keV, and a hydrogen column density of Nu ~ (6 ± 2) x 10 24 cm -2 . 
The non-thermal component can be fit with a slightly steeper photon index of r 2 = 2.0j^ 5- Adding 
a reflection component to the model (e.g. pexrav) does not lead to an improved fit. 

For the complete SED (4.5 keV - 2.5 GeV), we apply the model wabs * (gausslines + 
powerlaw + mekal ) +( wabs * cutoff pi ) +powerlaw. The first expression describes the 
absorbed soft X-ray energy range with a complex of five iron and nickel lines (gausslines 
where every single line is modeled by a Gaussian function), a photon scattering component (power law) 
and thermal starburst component (mekal). We fixed the line energies of the fluorescence emis- 
sions (Fei Kai, Ka 2 , K/3; Fexxv Ka; Nil Ka) to their expected positions according to the 
redshift of NGC 4945. Their positions are consistent with the values found when fitting the XIS 
spectrum only (Table §). The second expression is fitting the transmitted photons with a heavily 
absorbed exponential cut-off power law (wabs * cuto f f pi) and has the same photon index as the 
soft power law component. The third expression describes the non-thermal shock component in the 
y-rays as explained above. We can reproduce the spectrum of the AGN in the X-rays, the starburst 
(at soft X-rays), and the non-thermal shock component in the y-rays, achieving a fit statistics of 
X 2 = 1.16 (98 d.o.f). The combined data fits are shown in Figure [I] and the fit results are sum- 
marised in Table ||, including the iron and nickel lines detectable in the Suzaku/XIS spectrum. We 
obtain a total absorption corrected flux of ^keV-toOGeV ~5 x 10 _10 ergcm _2 s _1 and a luminosity 
of L 2ke v-iooGeV ~ 9 x 10 41 ergs" 1 . 



4. Discussion 

The model we applied to the combined data allows a representation of the spectral energy 
distribution over the whole 4.5 — 10 6 keV energy range we studied. For the X-ray range, we applied 
the models for the AGN and starburst at the same time. We obtain following fitting results: a photon 
index of T = 1 .6, a cut-off energy of E cut ~ 150 keV, an absorber in the soft X-rays (surrounding star 
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Figure 1: Left panel: unfolded spectrum E 2 f(E), Right panel: Soft X-ray spectrum with Gaussian lines 

burst component) with A^ oft = 7 x 10 22 cm~ 2 and A^ ard ~ 6 x 10 24 cm~ 2 for the AGN core emission 
itself. For the starburst we obtain a temperature of kT = 3 keV. These results are consistent with 
previous analysis of the Suzaku data alone [|6p, who used exactly the same data set but an earlier 
software and calibration version. In the y-ray domain the detection significance of NGC 4945 is 
still rather low. With the data we obtained of the AGN, we can reproduce the results of Lenain et 
al. [Q] which was based on an earlier Fermi/LAT data set. 

It is difficult to constrain the model for the combined /-ray and hard X-ray data, because of 
the wide gap in the ~ 200 keV to 100 MeV range. We applied two different models to explain the 
y-ray emission. In the first scenario, we assume that the y-ray and hard X-ray data are correlated 
and can be modeled by one non-thermal model. Here, the y-ray and the X-rays have to originate 
from the same inverse Compton process. In order to reach the y-ray domain, in this case one 
would need to assume that NGC 4945 hosts a jet, which gives rise to this non-thermal inverse 
Compton component. In the second model, we assume a thermal inverse Compton component for 
the hard X-ray spectrum of the AGN and a second non-thermal component, represented by the 
simple power law. Because of the low significance and relatively small bandwidth of the detection 
of this component, we can only assume that the y-rays are caused by the star burst in this case, 
for example by shock processes. Similar emission has been seen by FermilLKF in the star forming 
galaxies NGC 253 and M 82. One can think of this as a superposition of many Eta carina like 
systems, where a massive stellar wind hits the interstellar medium and thus produces strong shocks 
which give rise to the y-ray emission. Alternatively, cosmic ray interactions in the interstellar 
medium of the host galaxy can be responsible [f7|]. 

We obtain the best fit results with the second scenario, i.e. two independent components. Thus, 
it is more likely that the high-energy spectrum consists of a combination of a thermal (inverse 
Compton) and non-thermal component. NGC 4945 hosts no blazar, there is no jet observed in 
the source, and there is no strong radio emission which would be an indication for one single 
non-thermal component responsible for the hard X-ray spectrum. Also the presence of the strong 
iron and nickel fluorescence lines confirms the thermal and typical Seyfert-type AGN emission. 
The photons of the AGN disc are most likely the main source for the ionization and fluorescence 
processes in the AGN surrounding matter. 
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For the spectrum in the 2 — 10 8 keV energy range, we obtain a total absorption corrected flux of 
^2keV-iooGeV ~5x 10~ 10 erg cm~ 2 s _1 and a luminosity of L2keV-iooGeV ~ 9 x 10 41 ergs -1 . This 
flux and luminosity are similar to the value in the energy range of 2 - 300 keV. Thus, the X-ray 
energy range remains a good proxy for the total high-energy emission of AGNs, even if we consider 
the more complex case with y-ray emission as observed in NGC 4945. 

Assuming that the high-energy emission is about half of the total bolometric emission, we 
obtain an Eddington ratio of A ~ 0.005 for NGC 4945. Seyfert 2 galaxies normally tend to have a 
lower average Eddington ratio than Seyfert 1 galaxies ((Asy2) = 0.02 compared to (As y i) = 0.06, 
[§]). Other authors even find lower Eddington ratios, e.g. X = 0.006 for 87 AGN at z < 1.25 as 
presented by Simmons et al. On the other hand, low mass black holes accrete on average 
at higher Eddington ratios [TO]. Peculiar Seyfert galaxies can have Eddington ratios close to the 
Eddington limit (X = 1), like MCG-5-23-16 JiH], or even display super-Eddington accretion with 



X > 1 [12]. In summary, it appears that the Seyfert core in NGC 4945 is operating at a rather 
"standard" or low Eddington rate. The complexity of the X-ray spectrum with the transmitted, 
reflected, and star-burst component, is mainly due to the high signal-to-noise one can obtain for a 
near-by Seyfert galaxy. In the case of NGC 4945 the y-ray emission is not linked to the central 
engine, thus Seyfert galaxies remain to be no y-ray emitters, unless one considers the star burst in 
the host galaxy. The pre-requisite for y-ray emission from the AGN remains to be the existence of 



a powerful jet, like in the case of the y-ray bright source Cen A [13]. 
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